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Magic state distillation (MSD) is a purification protocol that plays a central role in fault tolerant quantum 
computation. Repeated iteration of the steps of a MSD protocol, generates pure single non-stabilizer states, 
or magic states, from multiple copies of a mixed resource state using stabilizer operations only. Thus mixed 
resource states promote the stabilizer operations to full universality. Magic state distillation was introduced 
for qubit-based quantum computation, but little has been known concerning MSD in higher dimensional qudit- 
based computation. Here, we describe a general approach for studying MSD in higher dimensions. We use 
it to investigate the features of a qutrit MSD protocol based on the 5-qutrit stabilizer code. We show that this 
protocol distills non-stabilizer magic states, and identify two types of states, that are attractors of this iteration 
map. Finally, we show how these states may be converted, via stabilizer circuits alone, into a state suitable for 
state injected implementation of a non-Clifford phase gate, enabling non-Clifford unitary computation. 



I. INTRODUCTION 

Quantum computers hold the promise of solving certain 
computational tasks at an exponentially faster rate than is cur- 
rently believed to be possible with classical computers flj. 
One of the main obstacles making the task of building a quan- 
tum computer difficult is due to quantum de-coherence, where 
errors can, if not corrected, accumulate and spread rapidly. 
The theory of quantum fault-tolerance (FT) |2| provides a 
mean to protect the coherent quantum state and allow a reli- 
able quantum computation. This can be achieved provided the 
physical error rate is below a certain threshold value. The ex- 
act value of the threshold error depends on the fault-tolerance 
scheme adopted for the computation. 

All FT schemes have a limited set of operations allowing di- 
rect implementation on the encoded quantum information J3), 
and in most known schemes these are the stabilizer operations. 
Stabilizer operations consist of a family of unitary circuits 
known as the Clifford group, preparation of |0) state and mea- 
surement in the computation basis. Although these operations 
can produce highly entangled states, the Gottesman-Knill the- 
orem @] tells us that a computation consisting of the stabi- 
lizer operations alone can be efficiently simulated by a classi- 
cal computer. The stabilizer state operations are not quantum 
universal, in fact, Clifford circuit operations alone cannot even 
implement non-linear classical logic gates, such as the Toffoli 
gate. Nevertheless, in most quantum FT schemes, the stabi- 
lizer operations are those which can be most readily achieved 
fault tolerantly. It is natural then to ask, what additional re- 
sources are needed to promote the stabilizer operations to uni- 
versality? 

One answer to this question is given by the theory of magic 
state distillation [5 6|. Almost any single qubit gate is approx- 
imately universal for SU(2) [7], and similarly, if the stabilizer 
operations are augmented by a supply of almost any ancillary 
pure states, the states can be used as resource for the imple- 
mentation of non-Clifford gates via "state-injection" methods 
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|6). A supply of qubits in a particular state may suffice, pro- 
vided it is not a stabilizer state. 

To compute fault tolerantly, we require a scheme that toler- 
ates preparation noise for these ancillary resources. Hence, we 
are interested in which mixed quantum states provide a suit- 
able resource. One can straight away identify a set of states, 
which will not be useful in this regard. These are the single 
qubit states that can be prepared via the stabilizer operations, 
together with classical randomness. These states are the con- 
vex hull of the Pauli eigenstates, and we shall call such states 
the stabilizer states. In the Bloch picture they occupy an octa- 
hedron, as shown in Fig. [T] 

Despite these obstacles, Bravyi and Kitaev \5\ showed that 
some mixed non-stabilizer states can enable universal quan- 
tum computing. They proposed a process of "magic state dis- 
tillation" whereby suitable resource states are efficiently con- 
verted, using only stabilizer operations, into a smaller num- 
ber of purer non-stabilizer states, the so-called magic states. 
While their methods are extremely powerful, they only con- 
sidered the qubit case. For non-generic noise models, it is 
also possible to distil certain 3-qubit states that can be used 
to implement the Toffoli gate [8|. Here we tackle the prob- 
lem of distilling pure 3-dimensional, or qutrit, non-stabilizer 
states, and herein refer to all pure non-stabilizer states as 
magic states. 

Known magic state distillation schemes have an iterative 
structure EID, with each iterate having 3 steps: 

1. Initialization, prepare n-copies of the qudit resource 
state p r ; 

2. Projection, measure Pauli-observables that stabilize a 
<i-dimensional subspace and postselect on the +1 out- 
come; 

3. Decoding, perform a Clifford unitary that maps the d- 
dimensional subspace onto a single physical qudit p' r . 

When successful, the output state p' r is used as one of the 
inputs on the next level of iteration. All known magic state 
distillation protocols achieve higher purifies via iteration, al- 
though non-iterative protocols are an interesting possibility; 
for example, compare with the hashing protocol iflOl [TP and 
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FIG. 1 : The Bloch sphere with the stabilizer octahedron and 
the magic states. There are 12 (blue) //-type magic states and 
8 (red) T-type magic states. 



quantum polar codes [ 12] used in analogous context of entan- 
glement distillation. If there exists a protocol that iteratively, 
or by other means, converts a state p r into a magic state of 
arbitrary high purity, then p r is said to be resource state. 

Magic state distillation for qubits has an elegant geomet- 
rical visualisation in terms of the Bloch sphere. In this rep- 
resentation, the stabilizer states form an octahedron, whose 
vertices are the Pauli eigenstates. Clifford Group unitaries co- 
incide with symmetries of this octahedron. The magic states 
distilled by the protocol of [6| correspond to the states invari- 
ant under certain rotational symmetries. The first of these 
is the set of 180° rotations around the edges of the octahe- 
dron, which contains the Hadamard gate, and we shall call 
Hadamard-type, or H-type rotations. The second type are 
120° rotations around the faces of the octahedron, known as 
T-type rotations. There has been a considerable amount 
of work to improve the original qubit schemes in [6] and to 
modify their noise model. Reichardt [ 13 1 showed how all the 
states above the edges of the octahedron can be distilled by 
Steane's 7 qubit code. In contrast, two interesting no-go theo- 
rems for qubit and odd-dimensional systems show that not all 
non-stabilizer states are useful resources. Recently, Veitch et. 
al. JT4 1 showed that for odd-dimensional systems there exist 
bound states that cannot contribute any enhancement to the 
computational power of the stabilizer operations. The qubit 
problem is more subtle. Campbell and Browne |[T5l showed 
that for any iterative protocol there will always exist undis- 
tillable qubit states above the faces of the octahedron. How- 
ever, non-iterative qubit protocols are still poorly understood. 
Furthermore, Campbell |[T6l introduced an activation protocol 
that can activate qubit states from above the octahedron face 
to the distillable regions. The magic state distillation protocol 
allows us to upper bound the error tolerance threshold ifTTl 
and is at the heart of state-of-the-art fault tolerance schemes 

ma. 

Surprisingly, very little study has been made of whether 
magic state distillation can be generalized to higher dimen- 
sions. There are a number of reasons why such a study would 
be important. In particular, topological quantum systems have 
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FIG. 2: An overview of how different protocols in this paper 
are related. We consider distillation of two different classes of 
states using the 5-qutrit code. The output of either can be con- 
verted by 2 subsequent sub-protocols, called parity-checking 
and equatorialization, to produce different magic states. The 
output of equatorialization can then be exploited to perform 
a non-Clifford unitary. The relevant section is noted in the 
top-right corner of each protocol box. 



anyons with braiding statistics and a dimensionality that is 
fixed by the physical system. Research into understanding 
these systems is ongoing [20, 24], but it is entirely plausible 
that the most promising systems provide braiding statistics 
corresponding to Clifford operations of dimensions greater 
than 2. Furthermore, the original magic state schemes have 
some surprising features (hence the name "magic states"), and 
it is unclear how much these depend on special features of the 
qubit Clifford group, or whether they are generic. There are 
therefore many open questions. Can the region of states which 
converge under distillation be identified? A study of higher 
dimensional MSD may give new insights into the structural 
differences between the Clifford group in two-dimensions, 
which plays a central role in quantum computation theory, and 
in higher-dimensions, where it has been intensively studied in 
the context of SIC-POVMS ED, MUBs J22]|2l and DWF 
(discrete Wigner functions) |24l . 

Recently, expanding on their work in l25ll . van Dam 
and Howard ll26ll have studied noise thresholds with d- 
dimensional quantum systems. They have calculated the set of 
robust qudit states that are most resilient to depolarizing noise 
and found the degree of noise needed to map such states to the 
set of stabilizer states increases with dimension scaling with 
d/(d + 1). Thus, the higher dimensional states have the po- 
tential to offer higher magic state distillation thresholds. Until 
now, however, no magic state distillation schemes for d > 2 
had been presented. 

In this paper, we present the first generalization of magic 
state distillation to higher dimensional systems. We demon- 
strate that magic state distillation can be achieved in a qutrit 
system, and find both similarities and differences with previ- 
ously known MSD protocols for qubits. We have chosen to 
focus on the three-dimensional qutrit space for this first study 
since it has the benefit of prime dimension, and computational 
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tractability. Nevertheless, many of the features which we un- 
cover are likely to be generic. 

The main protocol we consider is a qutrit generalization of 
the 5-qubit code covered in Sec. Ill and App. IB] Here we 



find magic state distillation occurring for 2 distinct families 
of states. The first family contains a pair of eigenstates of 
the qutrit Hadamard operator. Under depolarizing noise, they 
are distilled up to an error threshold of 23.3%. This fam- 
ily is generic in that all quantum states can be mapped into 
the Hadamard plane by random application of the Hadamard 
unitaries, a process known as Hadamard-twirling. The sec- 
ond family contain four eigenstates of the qutrit Hadamard- 
squared operator, but lies within a degenerate eigenspace of 
this operator, and so is not uniquely defined by it. Under de- 
polarizing noise, they are distilled up to an error threshold of 
34.5%. 

Not all magic states, such as those output by the 5-qutrit 
code, can directly perform injection of a non-Clifford gate. 
However certain special magic states, which we call phase 
states, do have the capacity to implement non-Clifford uni- 
taries. We show that, by executing additional protocols, the 
outputs from the 5-qutrit code can be converted into the desir- 
able phase-states. We overview the whole process in Fig. d2l, 
which shows that the promotion of the Clifford group can be 
achieved by following 5-qutrit distillation by parity-checking 
and equatorialization, which we describe in Sec. [V] 



II. DEFINITIONS AND NOTATION 

A quantum dit, or a qudit, is a <i-level quantum system 
where d > 2. The corresponding quantum state p can be rep- 
resented by a positive semidefinite d x d matrix of unit trace. 
In the language of linear algebra, such a matrix can be decom- 
posed as a linear sum of a basis set. For our purposes we have 
adopted the Weyl basis (also known as the Heisenberg-Weyl 
basis). This is because the Weyl basis set is a natural general- 
ization of the conventional qubit Pauli operators which makes 
them very convenient in studying MSD. For the remaining of 
this paper we will assume that d is a prime dimension. 



A. Higher Dimensions and Stabilizer Codes 

The rf-dimensional single-qudit X and Z Pauli operators 
are defined as 11271 : 

X = f \(j + 1) mod d)(j\, Z= £ J \j) 01 , (1) 



J=0 



j=0 



where u) = e 27Tl / d is the d-th root of unity. From this definition 
we see that X and Z are traceless non-Hermitian unitaries. 
They obey the commutation relation XZ = uj ZX and are 
cyclic in the power of d (i.e X d = Z d = I). We define a 
slightly different form of the single-qudit Pauli operators as 
follows: 



where c = (l-d)/2. The main reason for choosing this defini- 
tion and the significance of the extra phase oj c ^ k will become 
apparent in section II B But we shall first outline some of the 



properties of the Pauli operators based on this definition. For a 
single qudit, the composition of two Pauli operators can easily 
be verified to be: 

~ ~ , j'k—c(jk'+j k) „ /o\ 

cr^kCTj'M' = u °0+j"'),(fc+fc')- <- 3 ) 

For the case of a composite system of n-qudit we use the 
symplectic notation to represents the n-fold tensor products 
of Pauli operators: 

°hM ®<?hM ®' , -®°'3 n ^ = c r Ji:)2 ... i „,fc 1 fc 2 ...fc„ 

= &j,k, (4) 

where j and fe are vectors in ZJJ. The Pauli operators satisfy 
a generalized commutation relation: 



(5) 



where k.j - j.k' is the symplectic inner product. Based on 
the above definitions, the n-qudit Pauli group V n is: 



(6) 



All element of V n have eigenvalues of the form uj m for m € 

Z d . 

We will now briefly the stabilizer formalism and the basic 
properties of stabilizer codes [28 1. A stabilizer state \ip) is 
a simultaneous +1 eigenvector of an Abelian Pauli subgroup 
S n e V n called the stabilizer group. The subg roup S n is gen- 
erated by n independent and mutually commuting generators 
{51,52, ... , g n -k) with |<S„| = d n and -ui i S n . The stabilizer 
group S n forms a code, called the stabilizer code, with the 
stabilizer states being the codewords in the code-space. Stabi- 
lizer codes encodes k logical qudits (d k -dimensional Hilbert 
space Hk) into a larger Hilbert space H n of n physical qudits. 
We denote such code by [[n, k, 6]]d, where S is the distance 
of the code and subscript d is the dimension of the Hilbert 
space. Finally, we represent the logical operators on the code 
subspace by er^ fc , such that a^ k e V n \S n and commutes with 
all elements of V n . 



B. Qudit Space and Pauli Group Orbits 

When we discuss the general structure of a magic state dis- 
tillation in Sec. |V A| we will see that it is most convenient to 
use cr^fc as the basis set to represent a qudit state p. But recall 
that Oj^k is a non-Hermitian unitary operator. To guarantee 
the Hermiticity of p we must therefore impose the Hermitic- 
ity condition p = p' . We begin by expressing p as follows: 



p(a) 



(7) 



(2) 



where the summation is over all pair elements of l? d and we 
have assumed that a .o = 1- We will use a and the con- 
ventional 1 interchangeably. We refer to a j.k as the Bloch 
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components, as a generalization of the qubit convention, and 
a. as the Bloch vector which has the Bloch components as its 
elements. Observe that for d = 2, the set of cr^j. is Hermitian 
and the Bloch components will be real, but in the general qu- 
dit case the Bloch components are complex, constrained by a 
Hermiticity relation in order for p = p' to hold. To work out 
the Bloch components relation we start by explicitly writing 
P = P 1 ": 



E a 3,k< 7 j,k 



E 



(8) 



The importance of the extra phase factor of oj c ^ k in our defi- 
nition in eq. (|2|l is to ensure that k = C j.-fe, which can be 
easily verifier Using this fact, and after relabelling, eq. <|8j 
reduces to: 



-3-k- 



(9) 



A direct implication of eq. (|9]l is that only half the Bloch com- 
ponents, or (d 2 - l)/2, are independent, as the other half are 
simply the complex conjugates. Hence, only half the Bloch 
components are needed to define the density operator. Of 
course, the independent Bloch components are complex and 
we still have (d 2 - 1) real parameters defining the density op- 
erator. A Bloch component aj,k can be evaluated using the 
following relation: 



tr(pa + 



3,k 



.)=\x((XT-j- k ). 



(10) 



Beyond the qubit case, it is not possible to visualise the en- 
tire state space with a geometrical picture similar to the Bloch 
sphere. However, there have been some attempts to study the 
geometry of the state space in higher dimensions [30|, and 
for a Bloch type representation for qutrits BTl [32). Having 
fixed normalization and Hermiticity , a mixed qutrit state is 
described by a complex vector a, e C^ 1 ^ 2 . This complex 
vector space has a natural inner product, (a, (3) = Y,j a j Pj, a 
norm |a| = \J (ex, a), and distance D(a,(3) = \a - f3\. These 
geometric concepts are related to the density matrix represen- 
tation via 



HpIpp) = (l + 2(a,/3))/d. 
Since all pure state satisfy tr(p^) < 1, this entails 



a" < 



d-1 



(11) 



(12) 



and all physical states are within a Bloch-like ball of radius 
y (d - l)/2 about the origin, with pure states on the surface 
of the Bloch ball. The qubit state space is of course a special 
case in which all the points on the surface of the Bloch sphere 
corresponds to positive physical states. The additional condi- 
tion required to ensure a positive pure state, as shown in 11331 . 
is tr(p 2 ) = tr(p 3 ) = 1. 

The geometry vectors a corresponding to physical states 
p a is quite intricate. However, within certain hyperplanes the 
substructure is very simple. Consider hyperplanes defined by 
a set of d positive orthonormal operators, {pp 3 }, such that 



entails 



Note that in the geometric picture, this 



</3,,/3 fe ) = -(<%*-!) 



(13) 



We consider the complex hyperplane spanned by such 
that a = Y,j bjflj. It follows that an operator p a is positive 
if and only if Y,j bj < 1 and bj = \bj\ for all j. Hence, the 
physical a lie within the convex polytope with f3j as vertices. 
We have d vertices all equally separated from each other and 
residing within a real d-1 dimensional hyperplane. Hence, 
the polytope has the structure of a standard simplex. For d = 3, 
and an appropriate plane, the physical states reside within an 
equilateral triangle. 

Finally, we shall discuss the orbits of the Pauli group V n 
when acting on a general qudit state p(a) with conjugation 
being the group action. The orbits of a general state p{a), 
denoted by Orb(p(a)), is defined as: 



k' 



Orb(p(a)) = {p(a') = a f k - p{a) a 1 ., 

(14) 

In our Bloch representation, we are using Pauli group ele- 
ments as a basis set for the states, thus conjugation by Pauli 
operators will not transform the basis elements, but will add a 
phase of the form io l for some I e Therefore, the overall 
effect of this conjugation is to add certain phases to the Bloch 
components. The exact form of the phases is exactly given by: 



P( a ') = E a 3,fe (7 i',fe' CT i,fe a '-i',-fe' 
= u jk ~i k aj.k&j.k- 



(15) 
(16) 



where a ., ., = (Ji* _y, the commutation and composition re- 
lations where used in the last step. 



C. Qutrits 

So far, the discussion has been for all prime dimensions, 
but in the remainder of this paper we will discuss the qutrit 
case only. Therefore, we shall outline some of above results 
explicitly for the d = 3 case. Our definition of the qutrit Pauli 
basis set in eq. (teji is <j Jtk = lj'^X 3 Z k , where u = e 27 ™/ 3 and 
c = -1. The explicit qutrit p{a) state is: 

P( a ) = g^tho + ai.ooxo + a* 1Q a 2 ,o + ao,i<ra,i 

+ ao,i cr °> 2 + "M^M + "l.i ^ + "l, 2^1,2 + ai,2 CT 2,i)- 

(17) 

As we can see, completely specifying a qutrit state would only 
require 4 complex independent parameters. In terms of the 
Bloch components, the purity condition tr(p 2 ) = 1 for a gen- 
eral qutrit state can be shown to be \ot\ < 1. 

The Pauli group orbits for a single qutrit state can be eval- 
uated using eq. (16i. We are interested in knowing how 
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O"0,0 
ff±l,0 
ffO,±l 

ffil.Tl 



p(ai,o,ao,i,ai,i,ai,2) 

p(a; ±1 ai i o,^ Tl ao,i,ai,i! ajTla i.2) 
p(a; Tl Qi i o,^ Tl ao,i,ai,i w±1 ! a i.2) 



TABLE I: The qutrit orbital Bloch phases. 



the phases of the four independent Bloch components change 
when an element from the 9 qutrit aj f. operators is conjugated 
with the general qutrit state. The result is summarised in ta- 
ble [I] We refer to these phases as the orbital Bloch phases. 
These represent the phases which generate the set of states 
Pauli equivalent to any state. The magic states that we will 
find are unique up to a Bloch orbital phase. In other words, in- 
serting one of the phases from the set in table[j]into the Bloch 
components of the magic states would also give a valid magic 
state with the same distillation properties. 



D. Stabilizer Dynamics and the Clifford Group 

Indeed, one the most important advantages of the stabilizer 
formalism is the simplicity it provides when studying the dy- 
namics of the stabilizer states. Instead of studying the action 
of a unitary U on a n-qudit stabilizer state \ijj s ) (which would 
require the complete description of the map on d n parame- 
ters), the problem can be reduced to studying the evolution of 
the stabilizer operators by U in the Heisenberg picture (which 
is linear in n). An important class of unitaries are those that 
map the elements of the stabilizer group back to the stabi- 
lizer group under conjugation. These operation are called the 
Clifford unitaries. More formally, the Clifford group is the 
normalizer of the stabilizer group, defined as: 



C n = {C\CPi& = Pi V P it P, 6 V n ). 



(18) 



For any prime dimension, the Clifford group has been 
shown to be generated by three gates [27|[34]. These are the 
Hadamard gate H, the Phase gate S and the Controlled-NOT 
gate Apr). 
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V« keZ 



\k). 



(19) 



Under conjugation, the Hadamard gate transforms Pauli oper- 
ators as follows: 



(20) 



For completeness, the e?-dimensional S and A.(X) gates are 
given by: 

S\j) = ^-^\j), (21) 
Apr)|j)lfc)=|j)|(j + AO modd). (22) 



In addition to unitary dynamics, we allow for so-called 
Pauli measurements. For qubit systems the Pauli group con- 
tains Hermitian operators where for qudit they are unitary but 
not Hermitian, with the exception of the identity. Rather, 
when we say we measure a non-trivial Pauli P, this is taken to 
mean that we perform a POVM measurement that the eigen- 
vectors of P as POVM elements. 



III. MAGIC STATE DISTILLATION 

Using the definitions and notations we have developed in 
the previous section, we will show how the three steps of a 
MSD protocol described in the introduction can be formulated 
to study the distillation properties of any stabilizer code of any 
prime dimension. 

Resource state preparation: The computational model 
considered when studying MSD consists of perfect stabilizer 
operations and the ability to prepare n identical copies of a 
noisy resource state p r . It should be noted that the resource 
state is noisy due to the imperfect preparation process. By re- 
peating the preparation procedure n times, the state pf n will 
be prepared. As an input to the MSD protocol we consider a 
general state: 
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E 



jn.fcl — kn^jl— j n ,kl...k 



(23) 



The aim is, by performing the remaining steps of the proto- 
col on the above general form, we will determine the map on 
the Bloch components of the initial general state p. Then, by 
searching the state space for different initial states, we can 
identify the resource states as those that when used as an in- 
put to the protocol the output state has a higher fidelity with 
respect to a pure non-stabilizer state, and ultimately distilling 
this non-stabilizer pure state. If the search is done systemati- 
cally, one can in principle identify the entire region of resource 
states p r . 

Stabilizers measurement and Decoding: The (n-k) sta- 
bilizer generators of a stabilizer code [[n, fc, o]]d are mea- 
sured successively along with postselecting on the +1 out- 
come of each measurement. That is, if one of the outcomes 
is uj k (for some non-zero k e Z<j) then the protocol is aborted, 
and the procedure is repeated with a fresh state p %n . Also, 
it is important to notice that the error correction code is not 
being used for the usual purpose of correcting errors since the 
syndrome measurements are performed on the product state 
p® n . If successful, the measurement of the stabilizers simply 
project the state to the code's subspace. The projector oper- 
ator describing this measurement procedure can be put into a 
convenient form to us as follows: 



n = d fc 



k— r. 



mi 7rt2 

9i 92 



■9 n -V 



(24) 



After the measurements, the n copies of p will be projected 
into the code's subspace, and the following map will be per- 
formed: 

^--^y^. (25) 



tr(p®nn) 
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The state is decoded via a Clifford operator j9]- In a Heisen- 
berg picture, the decoding operation maps logical operators on 
the code-space, to unencoded operators acting on a single qu- 
dit. The output Bloch components after decoding q™! there- 
fore corresponding to the components of the logical operators 
prior to decoding a^ k . After one round of the distillation, 
these can be evaluated as follows: 



out 



trtiy^ny^) 

tr(p®«n) 



(26) 



The resultants expressions for the output Bloch components 
will be multi-variable complex polynomials of order n. For 
the qutrit codes we consider here, we have not found analytic 
solutions for the fixed points of the map. However, the prob- 
lem is tractable by using numerical methods to study the dis- 
tillation behaviours and the fixed points to a high accuracy. 



IV. THE 5-QUTRIT CODE 

Using the generalized formulation of MSD in the previous 
section we have studied the distillation properties of the five 
qutrit code. The stabilizer generators of the general five qudit 
code [[5, l,3]]d takes the same form in all dimensions. It 
is usually presented in terms of the conventional generalized 
Pauli operators of eq. ([TJ, as shown in table [II] 
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TABLE II: The stabilizer generators and logical operators of 
the five qudit code [[5, 1, 3]]^ 



Based on these stabilizers we have studied the distillation 
map of eq. ( B5 1 for the four Bloch components of a general 
input qutrit state. The exact distillation calculation and the 
expressions of the distillation map are given in appendix [B] 

The decoding of a stabilizer code is not unique, but one of 
an equivalence class of unitaries, a coset of the Clifford group, 
which are all equally valid choices. The choice of decoding 
will affect the iterative distillation behaviour. The decoding 
specified by the logical operators in table [n] is the canonical 
one, though we found behaviour was simplified by following 
each iterate with the following additional Clifford unitary: 



/ 



R 



1 



\ 



1 2 
UJ id LJ 



(27) 



/ 



where this maps a Hermitian operator p a such that: 

{ai,o,ao,i,ai,ii a i,2} {ot\ 2 > a i,i> a o ,1j a i,o}- (28) 



Without this corrective Clifford one observes a cycling be- 
haviour throughout the distillation process, see App|A] 

We identify two qualitatively different families of states. 
Firstly those in the Hadamard plane, which satisfy HpH^ = p. 
Secondly, we investigate distillation of an interesting set of 
states outside the Hadamard plane. 



A. Hadamard-like Distillation 

In the qubit case, the eigenstates of the Hadamard gate are 
known to be magic states, distillable by the five qubit code 
[[5, 1,3]]2 l35l . Since the exact generalized form of the 
Hadamard gate is defined in eq. Q19) >, a good starting point 
would be to investigate whether the qutrit Hadamard eigen- 
states can be distilled by [[5,1,3]]3. We begin by outlin- 
ing some of the structural properties of the qutrit Hadamard 
eigenspace. In the matrix representation, the qutrit Hadamard 
is given by: 



H 



1 

71 



1 1 



1 \ 



1 z 

1 W W 



1 



(29) 



where w 27 ™/ 3 , and it has the eigenvalues (+1,-1, i). We la- 
bel the corresponding three eigenstates as (\H + ) , , \Hi)). 
The density operators of the eigenstates have the form: 



\H+i) (H+i\ = p(a,a,b,b) , 
|£f_i><ff_i| = p(6,6,a,o), 
\Hi) (Hi\ = p(c,c,c,c) , 



(30) 
(31) 
(32) 



where a = \ (l+ \/3), b = \ (l- \/3) and c = -|, are real 
parameters. This basis of pure states all lie on the hyperplane 
of operators of the p(x,x,y,y). Probabilistic mixtures of 
these states form an equilateral triangle, and as reviewed ear- 
lier, all points outside this triangle correspond to non-physical 
operators. Furthermore, any qutrit state can be projected onto 
the Hadamard plane by applying the following twirling oper- 
ation to each copy of the input state p: 



(33) 



This twirling operation maps the general Bloch components 
as follows: 



cti i and «! 2 « 



Re(ai !0 + a 0; i) 
2 

Re(a M + ai j2 ) 



(34) 
(35) 



As such, we are interested in studying the distillable regions 
within the Hadamard plane, i.e. the states corresponding to 
the points inside the triangle. 

In studying the distillable region in the Hadamard plane, it 
is also informative to chart out regions for which distillation 
is impossible by any protocol. Clearly, all stabilizer states are 
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undistillable, but the results of Veitch et al [14| prove undis- 
tillability of all qutrit states with a positive Wigner function. 
The numerically calculated positive-region is shown in Fig. [3] 
i) as the red region. 

Running the [[5, 1,3]] 3 distillation for all the remaining 
points as inputs states we have discovered that both the \H+) 
and \HJ) states are distillable, but that \H.i) is not an attractor. 
The distillable regions are enclosed by the blue dashed trian- 
gles. The states \H+) and are equally valuable as magic 
states because R\H ± ) oc \H T ), where R is a Clifford unitary 
(see Eq. 



27 1 



which perhaps also explains the symmetry in 
their distillation regions. 

The path of the distillation takes the form shown in Fig. [3]ii) 
where we have chosen the \H+) blue triangle as an example. 
The small black points are few examples of input states to the 
distillation, and the black lines represent the distillation paths 
toward the \H+) state. Notice how the distillation does not 
follow a straight line (along the magic axes) as in the qubit 
case. In fact, in analogy to the qubit case, the above plane is 
the Hadamard magic plane. The curved distillation path can 
be understood by studying how the noise of a resource state 
in the Hadamard plane is suppressed in different directions by 
the distillation. We start by considering a general state p& 
inside the triangle of the following form: 

Pa = (1 - ex - e 2 ) \H + ) {H + \ + Cl (H-\ + e 2 \Hi) {H,\ , 

(36) 

with ei + €2 < 1. For clarity, we write the Bloch components 
of the above state as p a (A a , £?a , Ca , Da ) ■ They can be cal- 
culated explicitly using eq. (10) : 

Aa =B a = i(l + v^-2v / 3ei-(3 + v / 3)e 2 ), 

C A =D a = -(l--v/3 + 2^ci-(3--v/3)e 3 ). (37) 

Since we know the general distillation map for any set of 
Bloch components (see App. IB), we can simply substitute 
the above expressions into eqs. ( |B11|B14) to evaluate the 
output Bloch components a out = (A°^Bl ut , C A "*, D"^). 
The output state is then p° A = p(a° wt ). We have numerically 
calculated the output e° ut and e™' to the first order term as fol- 
lows: 



(ei,e 2 ) 
2 (£1,^2) 



out 



(H. 



I out 

I Pa 



H-) 



* (0.38 + 0.09e 2 )ei, (38) 
0.77e 2 + 3.55e 2 ei- (39) 



The above expressions shows an asymmetric error suppres- 
sion in the t\ (along the \H+) - line) and e 2 (along the 
\H+)-\Hi) line) directions. The particular distillation paths of 
Fig. |3]ii) can be explained by observing the difference in the 
coefficients of e out (ei,0) and e out (0,e 2 ), where we see that 
in the distillation region of the \H+) state there is a stronger 
attraction toward the \HA state compared to the \H-) state. 

The above analysis shows that the performance of the 
[[5,1,3]]3 in distilling the qutrit Hadamard states is not as 
good as the qubit case where the 15 qubit code by [6| has an 
output error probability of e out « 35e 3 . This is to be expected 
given the similar performance of the five qubit code 1351 in 
distilling the H-type qubit magic states. 




FIG. 3: A representation of the Hadamard plane. The 
Hadamard eigenstates are the vertices of the equilateral tri- 
angle and lies on a circle of radius l/\/2. i) The red region is 
the numerically calculated stabilizer region. The dashed blue 
and the green triangles contains the states that are distillable 
by [[5, 1,3]]3 and [[7, 1,3]]3 codes, respectively, ii) and iii) 
shows the distillation paths for the \H+) state and the mixed 
states for the [[5, 1, 3]]3 and [[7, 1, 3]]3 codes, respectively. 
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z 


z 
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TABLE III: The stabilizer generators of the seven qudit code. 



The state is not distillable by [[5, 1,3]]3. In fact, 
this state belong to the family of states with maximally non- 
positive Wigner function [26|. As we can see this state is the 
furthest away from the stabilizer region in the Hadamard plane 
and to bring it to the stabilizer region would require a depo- 
larizing noise with an threshold error of 0.75% (i.e. d/(d+l) 
for d = 3). Whether such a state is distillable is still an open 
question. 

To improve the size of the distillation region we propose 
a qutrit version of the seven qubit code [[7, 1,3]] 2 distilla- 
tion proposed by lfT3ll . We start with the stabilizer genera- 
tors of [[7, 1, 3]] 2 code and by adding the (-1) power to the 
appropriate X and Z Pauli operators, we construct a gener- 
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FIG. 4: i) The log-log plot of the output error probability for 
very small e. ii) The success probability for the trivial syn- 
drome measurements for the case where the magic states \ip) 
and \H ± ) undergoing depolarizing noise. 



alized 7-qudit commuting stabilizer generators as shown in 



table III We repeated the distillation procedure for these set 
of generators for the case d = 3 (exact calculations are omitted 
here) and we have investigated its distillation capability in the 
Hadamard plane. We found that this code attracts towards the 
non-stabilizer segments of the line joining the \H+) and 
states with the distillation region enclosed by the green trian- 
gle in Fig. [3] In other words, the 7-qutrit code distils, not 
pure, but mixed states. Regardless, the protocol may be use- 
ful for bringing states into the region distillable by the 5-qubit 
code. The distillation path for the \H + ) state is shown in Fig. 
|3]iii). This code increase the distillation region as shown by 
the solid blue curve in Fig. [3] For example, a state between 
the solid blue line and the dashed blue triangle is first distilled 
by the seven qutrit code to a state within the dashed blue trian- 
gle, after which the [[5, 1, 3]]3 code is used to distil the \H ± ) 
states. 



B. Hadamard-squared subspace 

In this section, we introduce a second class of magic states 
distilled by the [[5, 1, 3]]3 code. The state in question is an 
eigenstate of the H 2 operator, but lies within a degenerate 
eigenspace for this operator, and so is not uniquely defined 
by it. The magic state considered here has the form: 



|v>> = o|0) + 6|l) + 6|2), 



where up-to 4 decimal places 



a 
b 



-0.1203 
0.7017. 



0.0272i, 



(41) 
(42) 



The equality of the |1) and |2) components follows from the 
H 2 symmetry. In the Bloch representation, the a. vector is, 



a = {0.3236, -0.4772, 0.5438, 0.6098}, 



(43) 



with the feature that all components are real being again re- 
lated to the H 2 symmetry. Our numerical analysis shows that 
these states are clearly attractor fixed points of the distillation 
protocol, by we do not have a closed form analytic expressions 
for them. As a consequence, we will not be able to determine 
analytically how the error is suppressed as we did in the pre- 
vious section for the Hadamard magic states. Nevertheless, 
we can still gain a numerical indication of how the error of 
\(p) states is suppressed. Lets start with a \(p) state undergoing 
depolarizing noise: 



Pde P = (1 - e) \<p) {<p\ + el/3. 



(44) 



For a sufficiently small e the distilled state will also be of the 
above form (i.e on the depolarizing axis). We can then calcu- 
late the output error probability as follows: 



1 - (<fi\ Pdep\<P) 



(45) 



In general, we expect that e out a e n for very small e. Therefore, 
the power n can be evaluated as the gradient of a log - log plot 
of e versus e, as shown in Fig. 0i), and we found that n >» 1. 

For completeness we include the success probability p s of 
the syndrome measurement. A successful syndrome measure- 
ment, where all outputs of the stabilizer measurements is +1, 



is described by the projector II given in eq. (B15 1. Hence, the 
probability of this measurement is simply: 



Ps 



tr(np® 5 n) =trO® 5 n), 



(46) 



which is given in eq. ( |B16| l for all sets of Bloch components. 
We have computed p s for both \ip) and \H±) undergoing de- 
polarizing noise as an input states to the distillation. A plot of 
p s is given in Fig. |4]ii). 

V. PROMOTING THE CLIFFORD GROUP 

An interesting practical problem is how to use the magic 
states, \H+) and \ip), or their Clifford equivalent states, to per- 
form injection of a non-Clifford gate. While their utility is not 
immediately apparent in their current form, additional sub- 
protocols can be used to prepare the so-called /?/iase-states that 
are useful for gate injection. The phase-states are so-called be- 
cause they only hold phase information, having the following 
form: 



\^) = \0) + e w \l) + e^\2). 



(47) 



In Sec. V C we show how to use these states for gate injection. 



(40) However, first we describe, in Sec V A the parity-checking 
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protocol, which is used to convert both \H+) and \<p) into the 
plus-state \^ + ) = |0) + These plus-states are then input 



to the equatorialization procedure, in Sec. V B which finally 
outputs a desired phase state. For an overview of how these 
protocols fit together, the reader may refer back to Fig. pi. 



A. The parity-checker protocol 

Here we introduce a simple qutrit distillation protocol that 
is very efficient against a specific type of noise, but vulner- 
able against another type of noise. However, both \H+) and 
\(p) have zero overlap with the "bad" noise term and so the 
protocol can be efficiently used to convert these states into a 
plus-state. Before beginning the iterative protocol some ma- 
nipulation of the input states \H + ) and \tp) is required: 

1. Preparation 1, randomly choose from the set of unitaries 
{1,H 2 } and apply; 

2. Preparation 2, apply X ' ; 

3. Preparation 3, randomly choose from the set of unitaries 
{1, S, S 2 } and apply; 

where S= |0)(0| + |1)(1| + w|2)(2|. For eigenstates of H 2 the 
first step is not technically required, but it is listed to increase 
the generality of the protocol. In particular this preparation 
procedure maps all quantum states to 

p(So,tj) = (l-?7o-5o)|* + )(* + |+5o|^)(*1+?7o|2)(2|, (48) 

where |* ± ) = (|0) ± |l))/\/2. Note that |*") is Clifford equiv- 
alent to the Hadamard eigenstate \H.j). For imperfect \H+ ) and 
\cp) states, with depolarizing noise e, we have 



'7o = e(|V>» + £o/3, 
= £o/3, 



(49) 
(50) 



where for the two magic states of interest c(\H + )) = 0.2113 
and c(\(p)) = 0.0152. The parity-checker protocol will ex- 
ponentially suppress the value of whereas 5q will linearly 
increase. However, this is not problematic as So can be made 
arbitrarily small via distillation by the 5-qutrit protocol. 

The iterative parity-checker is now fairly simple, on the (n+ 
l)th round we have 

1. Take two copies of p(5 n , rj n ); 

2. Measure the observable ZiZ^ and postselect on +1; 

3. Decode the state such that \j, j) -> \j); 

4. Use the output state p(S n+ i, r/ n+ i) as an input in the 
next iterate. 

It is straightforward to verify the iterative relations are 

Vn+l = Vn/Pn, (51) 
Sn+1 = S n (l-i] n -S n )/p„, (52) 



where p n is the success probability 

p„=(l + r?„(3r ?n -2))/2. 



(53) 



When the small noise component is zero, so Sq = 0, and the 
large noise is not too large, r] < 1/3, then r\ n vanishes ex- 
ponentially quickly such that rj n ~ (2r)o) n . Allowing for 
non-zero Sq, the protocol can be iterated for approximately 
n ~ log((5o) rounds before the 8 noise becomes problematic. 

Let us consider a concrete example. If we have a 1^) magic 
state with depolarization noise e = 10~ 8 , this is first prepared 
into a noisy plus-state with r/o ~ 0.0152 and 5q = 10~ 8 /3. The 
total noise, r\ n + 5 n , will decease for the first 3 rounds of parity 
checking. After the 4th round we have a plus-state with a total 
error of only 2.707 x 1CT 8 . This illustrates that high-fidelity 
plus-states can be prepared from high fidelity \ip) states in a 
small number of rounds. 



B. Equatorialization 

Here we describe a simple magic state protocol that con- 
verts the plus-states to the phase states. The phase states lie on 
a generalization of the qubit Bloch sphere equator, hence the 
term Equatorialization. This protocol is probabilistic but not 
iterative. We take two highly purified copies of a plus-state, 
We measure a 2-qutrit stabilizer operator and postselect 
such that we project onto the subspace spanned by: 



|0l) 

|l i ) = X 1 X 2 |0 i ) 
\2 L ) = X 2 X 2 2 \0 L ) 



(|0,0)+ W |1,2) + W 2 |2,1)), 
(|1,1)+^|2,0) + W 2 |0,2)), 
(|2,2) +W |0,1) + W 2 |1,0», 



and then decode onto a single qutrit. When successful this 
produces the following transformation: 

|* + ) 82 ^(|0) + |l) + (w + w 2 )|2))A/3. (54) 

Noticing that cj + cj 2 = -1, we find that the output is a phase 
state 

l*o,») = (|0> + |1} - |2»/V3. (55) 
Finally, we have a magic state of the desired form. 

C. Gate injection 

Let us being by considering a general phase state \$>g^}. 
Given such a magic state, and a second qutrit state in any state, 
we perform gate injection by measuring Z\Z\. Given 
measurement outcome uj k , we perform the decoding \x, y) -*■ 
\y + k). The uj k outcome effects a unitary, Uk,e,<p, on that 
is diagonal in the computational basis with eigenvalues 



a i0 p i( P\ 



U ,oj = (l,e™,e*<»), 



Ui, 
U 2 , 



(56) 



i<f> 
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Each unitary occurs with equal probability because the phase 
state contains no variation in amplitudes. Although random, 
a desired unitary can eventually be reached by repeated at- 
tempts, with each attempt corresponding to a step of a random 
walk on a manifold of phase gates with a toroidal topology. 
For the |$o,tt) state, the corresponding unitaries take a simple 
form. The closure of t/fc o w gives a group of order 4, upto a 
global phase, which is composed of Uk,o,n and the identity. 
For such a small group any desired unitary will be reached, 
with high probability, within a small number of attempts. 

For brevity, let us herein denote N = C/o,o,7r- Clearly, N 
is non-Clifford and so we label the group C+m as the sin- 
gle qutrit group generated by TV and the Clifford group. The 
Gottesman-Knill theorem no longer applies and so we know 
no method of simulating computations using these unitaries. 
However, does this provide a dense cover of SU(3)? We dis- 
cuss this question in App. [C] where we show that the group 
C+m is of infinite order. The size of the group, and that it 
contains basis changing gates of the Clifford group, make it 
highly plausible that it provides a dense cover of SU(3). How- 
ever, we presently do not have a complete proof. 



group. We have offered good, but inconclusive, evidence that 
this promoted Clifford group enables universal quantum com- 
puting. 

Finally, we list few open questions that merit future inves- 
tigation. Do analogous magic state distillation protocols work 
similarly in arbitrary dimensions (or maybe just prime dimen- 
sions)? Are the magic states listed here exhaustive (could 
there be other attractive fixed points)? Can we find a closed 
form for the H 2 eigenstate coefficients? Can we better un- 
derstand the relationship between a code, and the magic states 
associated with the code? Is the set of unitaries (Clifford + 
(1,1,-1)) approximately universal? Are their topological quan- 
tum computing models which realise all or part of the qudit 
Clifford group? 
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VI. SUMMARY 



Appendix A: Clifford equivalences and cycling behaviour 



In this paper we have demonstrated that magic state dis- 
tillation protocols do exist for higher dimensional systems. 
We have provided a generic formulation that can be used to 
study the distillation capabilities of any stabilizer code for any 
prime dimension. We have showed that the five qutrit code 
[[5, 1, 3]] 3 is capable of distilling the qutrit Hadamard eigen- 
states, implying that, in analogy to the qubit case, there exist 
a i/-type qutrit magic states. Under depolarizing noise the 
|if } ± are distillable up to a noise threshold of 23.3%, this is 
a higher threshold than the |(3 - \/6) » 18.3% achieved five 
qubit code 1 35 1 . We have introduced a seven qutrit code and 
showed how this code can improve the distillation resource 
region of [[5, 1,3]]3. Interestingly for the case of the seven 
qubit code, the distillation threshold of the Hadamard eigen- 
states is tight. And although the seven qutrit code does not 
distil the Hadamard eigenstates, the numerical results of sec- 
tion A, suggests that this tightness is retained in the qutrit case, 
where all the mixed states along the convex line of \H+) are 
distillable. An important difference with respect to the qubit 
case is the existence of one Hadamard eigenstate which is 
not an attractor to the distillation protocols studied, but which 
forms an unstable fixed point. This state is also the Hadamard 
eigenstate furthest from the stabilizer region. 

After searching the state space for distillable regions, 
we have found another set of distillable states outside the 
Hadamard plane. These states are not eigenstates of a non- 
degenerate Clifford operator, and whether they are distillable 
by any other stabilizer code or are unique to the five qutrit 
code, is an open question. Under depolarizing noise, these 
states are distillable up to a noise threshold of 34.4%. 

We have also shown how to convert the outputs of the 5- 
qubit code into other magic states. In turn these have been 
shown to provide a unitary gate that promotes the Clifford 



We have chosen a particular decoding for which the dis- 
tillation protocol has the simplest behaviour. If instead, the 
canonical decoding was used, without the addition of an R ro- 
tation, then purification would still occur, albeit between each 
iterate the output would cycle between different states. In the 
Hadamard plane, we would observe an oscillation between 
\H±), which is obvious since R\H+) = \H+). Whereas, for 
the \ip) state there is a more complex 4-cycle, illustrated in 
Fig. [5] such that for the distillation map for one iterate, £, per- 
forms £(\<fij)) = \<pj+i) and \<p) = \<px) = \<p 5 ). The4-cycling 
states are related by; 



\¥>2) 



&\<pi), 

H\<pi), 
R f R?\<Pi). 



(Al) 



However, by considering £'{p) = R£'{p)R} , this cycling be- 
haviour vanishes. Note also, that this cycling behaviour is not 
only seen for the pure states but for depolarized states, and so 
all of these states are distilled by the 5-qubit code. 



Appendix B: [[5, 1, 3]]3 Distillation 

Based on the generalized magic state protocol constructed 
in section 



V A we will compute the exact distillation map on 



the Bloch components after a single round of distillation of 
the five qutrit code. 

We start by expressing the stabilizer generators and the log- 
ical operators of the [[5, 1, 3]]3 in terms of the aj ; k operators 
as shown in table [TV] The input to the distillation protocol is 



5 identical copies of a qutrit state p. Using eq. ( 23 1, the input 



11 



m) « 




1/3 



(1 - e) |^> +e 1/3 



\<p) 



Stabilizer States 
e > 



Distillable States 
e < 34.5% 



FIG. 5: An illustrative picture of the cycling behaviour of \<p). 
a) Starting with a mixed state (light blue point) the protocol 
will increase the purity of the states while cycling between 
them and ultimately reaching the fix pure points (dark blue 
points) . b) The convex line between one of the cycling states 
\(p) and the all mixed state 1/3 with an accurate ratio of the 
noise threshold. 
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TABLE IV: The stabilizer generators and logical operators of 
the five qudit code expressed in the o-j.fc notation. 



state is expressed as: 



Using eq. ( |24] >, a successful measurement of the four stabi- 
lizer generators with outcome + 1 corresponds to the following 
projector: 



(B2) 



Substituting the stabilizer generators in table|lV]into the above 
expression and using the composition law in eq. ([3]), the pro- 



jector becomes: 

n 4?> 



(9l-«3),(~«4) ® CT (92-94),(9l) ® Cr (? 3 ),(-9l+<?2) 



(-91+94), (-92 + 93) ® ^(-92), (-93+94) 



))• (B3) 



We can simplify the notation of the above expression by using 
eq. (|4|), w hich get rid of the tensor product sign, as shown in 
eq. < |B15) , 

The distillation map in eq. 
form as follows: 



( 26 1 can be put into a simpler 



out 

*3,k 



trtn^-IT'-^)-'-) 
tr(p®«n) 



tr(p®"n) 



3,k 



(B4) 
(B5) 

, II]=0 and the 



where in the last step eq. ( 10 1, II = II' 
cyclic property of the trace were used. 

The remaining task is to substitute eqs. (Bl 1 and (B3 1 into 
eq. (B5 1 to calculate the distillation map on the Bloch compo- 
nents. 

Lets start by evaluating tr( ( o® 5 II). Recall that all the <Tj. k 
operators are traceless except for the identity operator 00,0 ■ 
Therefore, the only terms that will survive in tr(p 05 II) are 
the coefficients of the identity operator. We get the identity 
operator in p® 5 II when the <jj : k operators in p® 5 and the oyjj 
in II have the opposite subscripts (i.e. a j' ,k' = "o.o 

iff 

j = j' and fc = fc')- As a result, tr(p® 5 II) will be the sum of 
all the Bloch components that are the coefficient of the a^k 
operators such that the subscripts (j, k) are the negative of the 
subscripts in eq. ( |B3[ l. In fact, since the summation is over 
all the elements of the ring Z| it is possible to multiply all 
the subscripts by (-1) without changing the actual value of 
the summation. Hence, tr(p lS ' 5 n) can be compactly express 
as shown in eq. (B16 1. In a similar way, we can express the 



tr (pHa_j _ k ) for all the four logical operators. For example, 
in the case of evaluating the output Bloch component a°" * , 



eq. (B5 1 becomes: 



e i,o 



tr(^ 5 rLr-\o) 
tr(p® 5 n) : 



(B6) 



withtr(p® 5 ncr L 10 ) given in eq. dB17k. 

We have evaluated the expressions for the four output 
Bloch components. However, writing them out in terms 
of ctj^ notation is cumbersome. Therefore, for clarity, 
we will relabel the four qutrit Bloch components as fol- 
lows (ai.o, ctp.i, ai.i, 0:1.2 ) = (A,B,C,D). For example, 
tr(p® 5 n) is given in eq. (BI81, where the subscript r rep- 
resent the number of the distillation rounds with r = cor- 
responding to the initial input state. Furthermore, it can be 
shown that the resultant expressions for the four output Bloch 
components can compactly be expressed in terms of a single 



function. This function is given in eq. (B19i, and based on 
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this function the distillation map can be expressed as: 

A r+1 =F(A,B,C,D), (B7) 

B r+1 =F(B*, A, D,C*), (B8) 

C r+1 =F(A*,C,B,D), (B9) 

D r+ i =T{B* ,D* ,A* ,C). (BIO) 

These four expression represents the complete distillation 
map, as the remaining four components are simply the com- 
plex conjugates of these expressions. However, the above ex- 
pressions do not incorporate the additional corrective Clifford 



(see Sec. IV I, and will result to the cycling behaviour in App. 
|A| We need to ensure that the map in eq. ( |28] i is applied after 
every iteration. This can easily be achieved in our formalism 
by the appropriate relabelling as follows: 



which is the corrected distillation map. We can see that in or- 
der to calculate the fixed points of this map analytically, one 
would have to solve the above simultaneous complex multi- 
variable polynomials of order 5. It is known from the famous 
Abel-Ruffini theorem that there is no algebraic solution for a 
general polynomial of order five or above. Therefore, the best 
way to discover the fixed points of the distillation is through 
numerical means. We started with initial states p(A, B, C, D) 
for a certain Bloch components and computed the above ex- 
pressions for a number of iterations, and observed whether 
there is a convergence toward a fixed point. 



A r+1 =F(D%C,B*,A), 
B r+1 =F(C*,D*,A,B), 
C r+1 =J^(D,B*,C,A), 
D r+l =T{C\A*,D,B*), 



(Bll) 
(B12) 
(B13) 
(B14) 



II oi S °'(9l-'/3)l(92-g4)2(93)3(-9l+g4)4(-<;2)5,(-<j4)l(9l)2(-gi+92)3(-'32+g3)4(-<;3+g4)5 ■ 

h"(p II) =— 2_, a (gi-<?3)l(92-«4)2(«3)3(-<?l+94)4(-«2)5,(-'34)l(gi)2(-9l+«2)3(-<?2+lj3)4(-«3 + 94)5- 
Ol ^4 



tr(p® 5 II< ) =^E« 



(9l-93)l (92-94)2 (93)3 (-91 +94)4(-92)s, (-94+1) l(9l + l)2(-9l+92 + l)3(-92 + 93 + l)4(-l?3 + 94 + l)5- 



tr(/9® 5 Il) =1 + 10(|A r | 2 + l^l 2 ) (\B r \ 2 + \C r \ 2 ) + h{B 2 r A* r C* 2 + D 2 A* r 2 B* r + 

D r (A 2 r D r C* + B 2 C 2 ) + B* 2 (A r C 2 + C* 2 D* r ) + D* 2 (A 2 r B r + C r A* 2 )). 

T(A,B,C,D) ={Bl + lOB r (D r A* + B* r ) (A r C* + C r D* r ) + 5(A r C 2 \A r \ 2 + D 2 (A r B* 2 + C r ) + 

Af (b; 2 d; + c;) + d; 2 (a 2 . + d t c* 2 ) + |c r | 4 s;))/tr( jO 05 n) . 



(B15) 
(B16) 
(B17) 

(B18) 
(B19) 



Appendix C: The size of the promoted Clifford group 

Here we show that C+m contains an infinite number of uni- 
taries. We do so by considering a particular unitary, 

K = H.N. H.N. H, (CI) 

which is an alternating sequence of Hadamards and our non- 
Clifford gate. The eigenvalues, e 47rA " of K can be analytically 
calculated by Mathematica, such that 



Since both Ai and A2 are irrational numbers, every power of K 
gives a different unitary, and so the single qutrit group, C+n, 
is infinite in order. The irrationality of Ai , and hence also A2, 
can be shown by expressing it as Ai = ArcCos(l/\/3)/7r and 
then using a theorem from ll36l . 



Ai = ArcTan(y2)/7r, (C2) 
A 2 = -ArcTan(\/2)/7r, (C3) 
A3 = 1/2. (C4) 
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